Disorder-induced coherent scattering in slow-light photonic crystal waveguides 
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We present light transmission measurements and frequency- delay reflectometry maps for GaAs 
photonic crystal membranes, which show the transition from propagation with a well defined group 
velocity to a regime completely dominated by disorder-induced coherent scattering. Employing a 
self- consistent optical scattering theory, with only statistical functions to describe the structural 
disorder, we obtain an excellent agreement with the experiments using no fitting parameters. Our 
experiments and theory together provide clear physical insight into naturally occurring light local- 
ization and multiple coherent-scattering phenomena in slow-light waveguides. 
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Introduction - Photonic crystals (PCs) are periodic 
dielectric structures that can strongly alter the propaga- 
tion of light due to the interaction of coherent reflections 
from the constituent periodic features. In particular, 
semiconductor-based planar PCs (e.g. Fig.[Tj) are of great 
interest due to their ease of fabrication using standard 
etching and lithography techniques. PC structures can 
exhibit interesting phenomena such as light trapping on 
sub- wavelength spatial dimensions in high-quality cavi- 
ties [H, or engineered waveguide band dispersions with 
a vanishing group velocity [H, 0, Both of these ef- 
fects give rise to novel regimes of light-matter interac- 
tion and have broad applications in nano- and quantum- 
technologies. 

An ongoing challenge, is the theoretical description of 
real fabricated samples, including the generally unknown 
role of unavoidable structural imperfections, collectively 
termed "disorder." Indeed, it is now well established that 
slow-light PC waveguides suffer from significant losses at- 
tributed to scattering at disordered surfaces and other 
device imperfections [HI, 0]. Rigorous modeling of this 
generally undesired extrinsic scattering phenomena is es- 
sential for understanding the underlying physics of mea- 
surements and for eventually producing functional de- 
vices. However, since the spatial scale of a disordered 
hole interface is only around 3nm or less (for ~200nm 
holes), and because of the complexity of the structural 
disorder [7] , the theoretical description of light scattering 
presents enormous challenges. 

Previously, Hughes et al. [sj introduced a PC wave- 
guide model for incoherent extrinsic scattering in the 
slow-light regime. The incoherent theory predicted the 
average scattering loss, using a first- and second- order 
Born approach for the scattering events, with the inco- 
herent averaging performed over a large set of nominally 
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FIG. 1: Left: Schematic of the planar PC waveguide. A sec- 
tion of the structure is cut away to show the anti-nodes in 
the y component of the Bloch mode electric field. Right: The 
computed corresponding band structure. The dispersion re- 
lations for the two lowest-order waveguide modes are shown 
and the lowest-order (blue, solid) mode is used for the calcu- 
lations. The grey region at the left represents the continuum 
of radiation modes above the light line. 

identical samples (an ensemble average). The "perturba- 
tive" backscattering loss was predicted to scale with the 
group velocity Vg as l/'U^ and dominate over scattering 
into radiation loss modes (1/'^^ scaling). Similar scalings 
were also implied by Povinelli et al. ^] and by Gerace and 
Andreani [lO|, though the scaling details depend on the 
nature of the propagating modes and how they sample 
the disorder regions. These approximate loss-scaling re- 
lations have been confirmed experimentally by a number 
of groups, e.g. but they may break down at low 

group velocities where multiple disorder-induced scatter- 
ing becomes significant. Recent experimental measure- 
ments have reported interesting features such as narrow- 
band resonances near the band edge [1, [l^ that are not 
explained at all by the incoherent and perturbative scat- 
tering theories (see Figure [2]) and simple Vg scaling rules. 
In this Letter, we introduce a non-perturbative the- 
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ory of coherent optical scattering. In contrast to pre- 
vious works, we include unlimited forward- and back- 
scattering events (through a coupled mode theory) which, 
combined with modeling individual, fully 3D^ disordered 
waveguides instead of ensemble averages (a statistical av- 
erage over many disordered waveguides), is essential to 
understand the new experimental observations. The the- 
oretical results are presented along side measurements on 
state-of-the-art GaAs PC waveguides which are probed 
with transmission measurements 0] and optical low- 
coherence reflectometry (OLCR). OLCR allows the mea- 
surement of the back-reflected signal as a function of 
propagation time inside the PC waveguide. The exper- 
imental setup is basically a Michelson interferometer il- 
luminated with a broadband, partially coherent source. 
The waveguide is placed in one arm and a translating 
mirror in the other. Properties such as group velocity 
and complex frequency-dependent reflectance [l5[ can be 
extracted from the interference pattern. We also present 
a powerful and recently developed analytic technique: 
time- frequency reflectance maps (TFRMs) [6], that can 
be used to visualize the frequency-dependent impulse re- 
sponse and reveal a number of interesting features such 
as disorder-induced scattering and facet reflections. 

Our formalism also provides fresh insights into the long 
standing question of light localization in PC waveguides. 
It was proposed by John [l^ and Anderson |l7| . that 
strong localization may be observable in PCs. Although 
Anderson localization was originally proposed for elec- 
trons propagating in a disordered atomic lattice, the phe- 
nomena is general to waves in a periodic medium and ap- 
plicable to PCs. Anderson localization occurs when the 
mean free path of a propagating Bloch mode is reduced 
to the order of the Bloch wave vector /c, totally disrupt- 
ing propagation; formally kl < 1 [Tsj. There have been a 
number of milestones towards observing Anderson local- 
ization in PCs including the observation of localization in 
a purely random powder [1^] , and localization transverse 



to the propagation direction [20|, |2l| . Recently, Topolan- 
cik et al. observed resonances at the band edge of an ar- 
tificially roughened PC waveguide and argued these were 
due to localization [l^. In contrast, here we report mea- 
surements on structures with only unavoidable disorder^ 
and provide deeper insight into localization phenomena 
and the correct criteria for strong localization. 

Theory .- In the spirit of a coupled mode approach, 
we introduce "slowly- vary ing envelopes," ?/^f[b](x) for the 
forward [backward] wave, to approximate the solution as 

E(r; CO) =So[ek{r) e'^^ ^x) + e^(r) e'^^^ ^^,{x)] + . . . , (1) 

where . includes contributions from lossy radiation 
modes, £q is an amplitude, e/c(r) is the periodic Bloch- 
mode electric field normalized by J^^jjdr£^(r) |e/e(r)p = 
1, and k is the wave vector which implicitly depends on 
the the angular frequency uo. 



We begin with an ideal structure (no disorder), de- 
scribed through the spatially-dependent dielectric con- 
stant £i(r), as schematically illustrated in Figure [H The 
structure supports a Bloch waveguide mode (see the solid 
blue curve in the band structure plot), which exists below 
the light line (and is thus lossless in the absence of dis- 
order), and which has a group velocity Vg = duo/dk that 
tends to zero at the band edge. If the electric field of the 
ideal mode E^(r; cj), is incident on a disordered waveguide 
described by £{y) = £^(r) + A£(r), the exact electric field 
has a well known analytical form: E(r;C(;) = 'Ei{r;cj) + 
/dr' G'(r,r';cj) • [E{r';u;) As{r% where G{r,r';u;) is 
the Green function for the ideal structure. The Green 
function is obtained from a polarization-dipole solution 
to Maxwell's electrodynamics equations. By writing 
E(r; u) as in Eq. ([T]) and by decomposing the Green func- 
tion as the superposition of an exact bound mode con- 
tribution and a radiation modes background summa- 
tion G rad(r, r'; cc;), we then derive a pair of coupled- mode 
equations for the evolution of the envelopes: 



d^pfjx) 
^ dx 

' dx 



i Cff (x) V^f {x) + i Cfb (x) e ^^^"^ V^b {x) 

+iCfr(x) V^f(x), (2) 

^Cbb(^) i^h{x) +zcbf(x) e^^^^ V^f(x) 

+iCbr(:r) V^b(^). (3) 



The coupling coefficients can be physically interpreted as 
Cff = Cbb driving scattering from a mode into itself, Cbf = 
c^^ driving scattering into the counter-propagating mode, 
and Cfr and Cbr driving scattering from the waveguide 
mode into radiation modes above the light line. We have 
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(4) 
(5) 



dydzdr' Ae{r)e*^ f^{r) ■ G rad(r, r'; cj) 
en,fc(rOAs(rO, n = f,b, (6) 



where ef^/e(r) = ^(r) = ej.(r) e*^^, and a is the period- 
icity of the PC. The most interesting frequency response 
of the system is dominated by scattering between wave- 
guide modes, through Vg{u;); while radiative scattering 
merely leaks energy from the waveguide and is a much 
smaller effect in the slow-light regime as shown in previ- 
ous studies. Importantly, this theory thoroughly extends 
previous formalisms, e.g. [sl], to account for i) multiple 
scattering events, ii) coherent scattering, and Hi) a di- 
minished Bloch mode amplitude. In addition, we solve 
the coupled mode equations numerically for disordered 
waveguides "instances," in a fully self-consistent way. 

Experimental Device .- Except where stated other- 
wise, the experimental device is a Wl.l PC waveguide 
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fabricated from GaAs as schematically shown in Fig- 
ure [TJ The width of the waveguide is l.l>/3a, the pitch 
is a = 400 nm, the thickness is 265 nm, the hole radius 
is R = 0.27 a, and the length is 250 /im. The fabrication 
quality of our GaAs devices is comparable to state-of- 
the-art silicon processes [22| . The samples were analyzed 
using the high resolution SEM and the image processing 
technique of Skorobogatiy et al. [7]. The disorder was 
found to be well described by small perturbations of the 
radius around the hole perimeter AR{(j)a)^ that follow the 
distribution {AR{(t)a)AR{(t)'^^)) = (j'^e^\^--^'o^'\^^^5^^oc' ^ 
where a indexes the holes, and (p^ is the angular coordi- 
nate of the point. The RMS roughness a, and correlation 
length /p, are estimated to be 3 and 40 nm respectively 
and these values are used in our calculations. 

Transmission spectra .- An experimental transmission 
spectrum is shown in the top plot of Figure [2] for a 
1.5 mm Wl waveguide of a different design. Approaching 
the band edge, the transmission rolls-off approximately 
with 1 /'^^ , however there are numerous sharp resonances 
where the transmission varies by orders of magnitude. 
Two theoretical models for this waveguide are shown in 
the bottom plot. The previous incoherent loss calculation 
[1] (dashed, red) captures the approximate l/'i;^ roll-off 
but does not explain the resonances. In contrast, the new 
coherent loss calculation presented in this paper (solid, 
blue) reproduces them since it accounts for multiple scat- 
tering events which are necessary to build up Fabry- 
Perot-like resonance between disorder sites. For reference 
the group index is show in the inset (smooth, blue). 
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FIG. 2: Top: Experimental transmission spectra for a 1.5 mm 
Wl waveguide showing resonances near the band edge. Bot- 
tom: Theoretical transmission spectra calculated using the 
incoherent (dashed, red) and coherent (solid, blue) scatter- 
ing theories. Inset: Ideal group index ng — c/vg (smooth, 
blue) compared with an estimate of the effective value due 
to disorder (red, rough). The effective group velocity is soft- 
ened due to disorder and does not diverge like the ideal value. 
The experimental spectrum is shifted by 1 THz to account for 
uncertainties with the actual fabricated slab thickness. 



While we obtain the general trends of the experiments, 
over more than three orders of magnitude (and without 
any fit parameters), we have not included the disorder- 
induced broadening of the slow-light-regime band struc- 
ture which alters the effective group velocity [23] and 
softens the roll-off. From a simple perspective, disorder 
adds or removes dielectric from an unperturbed (to first 
order) Bloch mode, causing a local frequency shift of the 
band structure. A propagating mode at fixed frequency 
thus has a locally- varying group velocity due to disorder 
shifting the waveguide band in frequency. From pertur- 
bative calculations with identical disorder statistics, we 
estimate that the group velocity will be noticeably al- 
tered from the ideal value for Vg ^ c/35 and will have 
a minimum of around c/80 at the band edge, as shown 
in the inset of Figure [2] (rough, red). For different struc- 
tures, this minimum Vg will vary. Although these findings 
are not important for the reflectance maps analyzed be- 
low, our calculations are broadly consistent with our own 
experiments, those reported by Engelen et al. [lH, and 
theoretical analysis of PC coupled-cavity structures 0. 

Time- Frequency Reflectance Maps.- TFRMs are in- 
tensity plots of the reflected signal as a function of time 
(horizontal axis) when the structure is excited with a 
narrow-band pulse centred at some frequency of interest 
(vertical axis). The map is generated using the complex 
reflectance of the waveguide which, for physical samples, 
is deduced from a single set of OLCR data or, for simu- 
lated structure, is calculated directly [6]. 

A simulated TERM is shown in Figure [3] for a wave- 
guide with perfectly transmissive facets. The dashed blue 
lines indicate the time the pulse is injected and the round- 
trip time. The solid magenta line indicates the ideal 
group index Ug = c/vg for comparison. Away from the 
band edge, at higher frequencies, the back reflections are 
small and are confined between the two time limits, in- 
dicating that only single scattering events are occurring. 
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FIG. 3: A TERM for a 250 /xm simulated disordered waveg- 
uide with perfectly transmissive facets showing the strength 
of the back reflection in dB. The left blue dashed line indicates 
when the pules was injected and the blue dashed curve is the 
expected round trip time in an ideal structure. The magenta 
line shows the group index Ug — cfvg on the top scale. 
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FIG. 4: Comparison of simulated (top) and experimental 
(bottom) TFRMs for a 250 /xm waveguide with a partially 
reflective front facet and a strongly reflective back facet. The 
simulation uses the same disorder configuration as in Figure [S] 

Approaching the band edge, for Vg < c/20, scattering be- 
comes significant with strong back reflections and multi- 
ple scattering events clearly evidenced by the hot spots 
and the continued reflections after the time for one round 
trip. This agrees well with Engelen et al. [l2[ who ob- 
served total disruption of propagation for Vg ^ c/30. 

Experimental samples have more complicated TFRMs 
due to reflections from the sample facets. Figure [H 
compares simulated (top) and experimental (bottom) 
TFRMs for the same device geometry. The multiple facet 
reflections [26*] are clearly visible and the time for a round 
trip lengthens as the group velocity decrease. At low 
group velocity, the pulse is washed out by strong mul- 
tiple scattering events, making transmission of signals 
difficult. There is an excellent correspondence between 
the measured and simulated maps. Since the simulations 
do not incorporate details of the experimental setup be- 
yond the sample, they tend to give richer features than 
the experimental maps. Nevertheless, there is an excel- 
lent correspondence between the measured and simulated 
maps^ and this agreement has been found for a number 
of different sample lengths and facet reflections. 

Localization .- Both our measured and simulated 
transmission spectra (Figure[T]) exhibits sharp resonances 
near the band edge, similar to those reported in [13] • 
These features can also be resolved in high resolution 
TFRMs. To rigorously address the question of whether 
these features are indicative of localization, a localization 
length / can be defined as l~^ = —{lnT)/L where T is 
the transmitted power and L is the sample length 
For the experimental structure at = 0.45 x 27r/a where 
Vg = c/4:b, the localization length is calculated to be 
~ 100 /im, far from the criteria for strong localization. 



Experimental imaging of light leaking from the waveg- 
uide plane confirms that these features are distributed 
over a large number of waveguide periods. Thus these 
features are better described as Fabry- Perot- like reso- 
nances between scattering sites and not localization, in 
agreement with the interpretation of Vlasov et al. [25[. 
This in no way alters the fact that backscattering in slow 
light modes leads to highly disordered propagation and 
low transmission, which has important implications for 
fundamental optical physics and slow light applications. 

Conclusion - By matching several carefully-designed 
light propagation experiments with a new non- 
perturbative, coherent scattering formalism, we have 
demonstrated and explained the crossover from nominal 
light propagation, with a well defined group velocity, into 
a naturally disordered regime that is dominated by co- 
herent scattering and light localization. Our simulated 
waveguides show excellent qualitative agreement with 
measurements on GaAs waveguides, as demonstrated by 
comparing the rich features of the TFRMs. We also de- 
termine the localization length near the band edge of a 
PC waveguide and show that while propagation is highly 
disordered, strong localization is not occurring. 
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